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Abstract 


Biochar management has been proposed as a possible tool to mitigate anthropogenic CO, emissions, and thus 
far its impacts in forested environments remain poorly understood. We conducted a large-scale, replicated field 
experiment using 0.05-ha plots in the boreal region in northern Sweden to evaluate how soil and vegetation 
properties and processes responded to biochar application and the disturbance associated with burying biochar 
in the soil. We employed a randomized block design, where biochar and soil mixing treatments were established 
in factorial combination (i.e., control, soil mixing only, biochar only, and biochar and soil mixing; n = 6 plots of 
each). After two growing seasons, we found that biochar application enhanced net soil N mineralization rates 
and soil NHj concentrations regardless of the soil mixing treatment, but had no impact on the availability of 
NO, , the majority of soil microbial community parameters, or soil respiration. Meanwhile, soil mixing 
enhanced soil NO, concentrations, but had negative impacts on net N mineralization rates and several soil 
microbial community variables. Many of the effects of soil mixing on soil nutrient and microbial community 
properties were less extreme when biochar was also added. Biochar addition had almost no effects on vegetation 
properties (except for a small reduction in species richness of the ground layer vegetation), while soil mixing 
caused significant reductions in graminoid and total ground layer vegetation cover, and enhanced seedling sur- 
vival rates of P. sylvestris, and seed germination rates for four tree species. Our results suggest that biochar 
application can serve as an effective tool to store soil C in boreal forests while enhancing NHj availability. They 
also suggest that biochar may serve as a useful complement to site preparation techniques that are frequently 
used in the boreal region, by enhancing soil fertility and reducing nutrient losses when soils are scarified during 
site preparation. 
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management has been proposed as an approach that 


I : 
a may potentially accomplish these goals simultaneously, 


Since the start of the industrial revolution, atmospheric 
CO, concentration have risen from 280 to above 
400 ppm, with some models predicting concentrations 
will reach 1000 ppm by the end of the century (IPCC, 
2013). These emissions have already led to increases in 
mean global temperatures and caused associated 
changes in precipitation in many regions (IPCC, 2013). 
Due to the current and potential impacts of these 
changes on society, there is substantial interest in devel- 
oping sustainable land management systems aimed at 
sequestering atmospheric carbon (C) while simulta- 
neously maintaining and optimizing other ecosystem 
services such as commodity production (Burney et al., 
2010; Palm etal., 2010; Post etal., 2012). Biochar 
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and as such has gained widespread attention among 
researchers and policy-makers throughout the world 
(Lehmann & Joseph, 2009). Biochar management 
involves intentionally charring biomass and burying it, 
with the purpose of directly storing C in the soil, while 
at the same time enhancing soil quality and plant pro- 
ductivity (Lehmann, 2007a,b; Sohi, 2012). It has been 
proposed that a significant proportion of anthropogenic 
CO, emission could be mitigated if biochar manage- 
ment were adopted globally (Lehmann, 2007b); how- 
ever, key data are lacking regarding its impacts on soils 
and plants in a variety of ecosystem types (Sohi, 2012; 
Jeffery et al., 2013; Cusack et al., 2014). 

The vast majority of biochar research has to date been 
conducted in tropical and temperate agro-ecosystems 
(Atkinson et al., 2010; Jeffery et al., 2011; Lehmann et al., 
2011; Clough et al., 2013). Studies from these systems 
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have shown that a wide variety of charring techniques 
can be used to stabilize a range of organic residues, by 
altering their C chemistry (e.g., structural components 
such as cellulose or lignin) into irregular or graphitic 
molecular structures (Steiner et al., 2007; Nelissen et al., 
2012; Clough et al., 2013). Many microbes lack efficient 
enzyme systems to mineralize charred materials, allow- 
ing some charred materials to remain stable over cen- 
tennial to millennial time scales (Glaser et al., 2001; 
Preston & Schmidt, 2006). While it’s very high C content 
(usually >70% C) and recalcitrance have the potential to 
contribute to long-term soil C storage, numerous studies 
have also shown that biochar can enhance soil fertility 
and plant growth. For instance, biochar often has high 
concentrations of available nutrients (e.g, NHj , PO; , 
Ca*’, and Mg*’) on its surfaces, which can have fertil- 
ization effects over short time-scales (Gundale & Delu- 
ca, 2006; Chan & Zhihong, 2009; Jeffery et al., 2011; 
Pluchon et al., 2014). Biochar has also been shown to 
enhance nutrient availability over longer time scales by 
enhancing nitrogen (N) mineralization or nitrification 
(DeLuca et al., 2006; Ameloot et al., 2015) as a result of 
enhanced microbial growth and activity (Zackrisson 
et al., 1996; Lehmann et al., 2011) and by reducing soil 
nutrient losses due to its high ion exchange capacity 
(Atkinson et al., 2010; Brewer et al., 2011). Numerous 
recent studies have shown that the positive effects of 
biochar on soil fertility can result in enhanced plant 
growth (Wardle et al., 1998; Lehmann, 2007a; Jeffery 
et al., 2011; Biederman & Harpole, 2013; Liu et al., 2013), 
thereby having an indirect positive effect on net ecosys- 
tem C uptake. 

A very limited amount of research on biochar has 
been conducted in forested ecosystems, including boreal 
forests that cover approximately 11% of the terrestrial 
land surface area and serve as one of the largest stores 
of terrestrial C (Sabine et al., 2004). A wide range of 
intensive forest management approaches are currently 
being investigated in forest ecosystems to either seques- 
ter C or produce biofuel to offset fossil fuel demand 
(e.g., stump harvesting, establishment of short rotation 
age exotic tree plantations, and biochar management) 
(Walmsley & Godbold, 2010; Lindholm et al., 2011; Sa- 
batti et al., 2014). One aspect of biochar management 
that remains poorly understood in forested systems is 
whether biochar can enhance soil fertility without stim- 
ulating loss of existing soil C. Several aspects of biochar 
management have the potential to cause unintended C 
losses, including physical disruption of soil during 
burying (Jones et al., 2010) and positive effects of bio- 
char on soil microbial growth and activity (Pietikainen 
et al., 2000; O'Neill et al., 2009; Liang et al., 2010). Previ- 
ous studies have suggested that the addition of charcoal 
to boreal soils can accelerate the decomposition of 


organic matter (Wardle et al., 2008); however, studies 
have yet to discern the effects of biochar itself vs. those 
of physical disturbance on soil C losses. An additional 
aspect of biochar management that remains poorly 
understood in forested ecosystems is its potential 
impact on ground layer vegetation composition, as well 
as the regeneration of forest tree species (Pluchon et al., 
2014). In contrast to most agricultural systems, a wide 
variety of plant species often co-occur in forested envi- 
ronments and these can exhibit considerable variation 
in their functional traits (Hobbie, 1996; Cornwell et al., 
2008) and relationships with local-scale soil fertility 
(Wardle et al., 2012). Thus, fertility changes associated 
with biochar application may impact ground layer spe- 
cies composition, which could in turn influence tree 
seedling establishment and growth rates. Such effects 
are relevant to understanding ecosystem C sequestra- 
tion because as forests develop trees eventually serve as 
the largest aboveground C pool. To our knowledge, the 
impact of biochar management on soil nutrient and C 
dynamics, as well as ground layer vegetation composi- 
tion and tree seedling establishment, has not been eval- 
uated at an operational scale in a forested environment 
(Robertson et al., 2012; Sohi, 2012; Cusack et al., 2014). 
We conducted a large-scale experiment in the boreal 
forest zone of northern Sweden to understand how bio- 
char management impacts on a wide range of soil and 
plant properties. We applied biochar and soil mixing to 
field plots in a factorial combination (i.e., control, bio- 
char only, mixing only, and biochar and mixing com- 
bined) so that we could discern the impact of biochar 
addition by itself (ie., applied on the soil surface) vs. 
the disturbance impact caused from incorporating it 
into the soil. We hypothesized that: 1) the addition of 
biochar will enhance nutrient availability and minerali- 
zation rates, microbial biomass, and soil respiration 
rates independently of the soil mixing treatment; 2) soil 
mixing alone also will enhance nutrient availability and 
mineralization rates, and soil respiration, but decrease 
total microbial biomass and fungal to bacteria ratios, 
and further, that soil mixing will enhance any effects of 
biochar on these properties through bringing the bio- 
char in closer contact with the soil; 3) enhanced short- 
term nutrient availability caused by the individual and 
synergistic effects of biochar, and mixing treatments 
will result in significant changes in ground layer plant 
species composition, with higher nutrient availability 
leading to greater plant cover and dominance of grami- 
noid species, which have high relative growth rates; 
and 4) greater ground layer vegetation cover associated 
with biochar and mixing treatments will correspond 
with reduced germination of tree seeds and survival 
rates of tree seedlings, due to competitive suppression. 
Testing these hypotheses in combination will provide 
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valuable insights into the impacts of biochar manage- 
ment on ecosystem functioning in forest ecosystems. 


Materials and methods 


Study site and experimental design 


The experiment was performed at the Aheden research area 
within the Svartberget Experimental Forest (64°14'N, 19°46’E, 
175 m above sea level) in northern Sweden. Soil at the site is a 
fine sandy Typic Haplocryod (FAO, Cambic Podzol) formed 
from silty glacial outwash sediment. The annual mean air tem- 
perature at the site is + 1.0 °C. Snow usually covers the ground 
from the end of October to late April. Mean annual precipita- 
tion is approximately 600 mm, of which half falls as rain and 
half as snow (Gundale et al., 2011a). Prior to the start of the 
experiment, the experimental site was covered with a closed 
tree canopy consisting of ~60-year-old Scots pine (Pinus sylves- 
tris). Prior to the experiment, the ground layer vegetation con- 
sisted primarily of ericaceous shrubs, mainly Vaccinium vitis- 
idaea and Calluna vulgaris, and of mosses and lichens (predomi- 
nantly Pleurozium schreberi, Dicranum sp., Cladina rangiferina, 
and Cladina arbuscula (Gundale et al., 2011a). 

In October, 2010, a 5-ha area within the experimental forest 
was clear felled using standard forestry practice, and tree boles 
and branches were removed from the site. The following sum- 
mer (i.e., July 2011), soil C stocks were estimated at the site by 
measuring the dry mass and C concentration on 24 homoge- 
nized soil cores, which were 20 cm in depth and included the 
organic horizon and mineral soil. These measurements 
revealed that the average bulk density across these soil depths 
was 1.00 g cm, and in total contained 38.8 Mg of C ha! 
(standard error = 2.82) prior to the establishment of treatments. 
In October, 2011, we established 24 0.05-ha (22.4 x 22.4 m) 
experimental plots, with a 5-meter buffer zone between each 
plot. The plots were arranged into six blocks, with each of four 
treatments (control, biochar addition only, soil mixing only, 
and biochar addition followed by soil mixing) randomly 
assigned to one of four plots within each block. Because 
stumps were an obstacle for effectively implementing the soil 
mixing treatment, we used an excavator which we positioned 
outside of each plot to vertically lift and remove all stumps 
from every plot. This removal caused a local disturbance at the 
site of each stump, but left the organic horizon intact in most 
of the plot. We note that stump harvesting is an active area of 
research in the Nordic region because of its considerable poten- 
tial as a source of biofuel (Walmsley & Godbold, 2010; Lind- 
holm et al., 2011). Following stump removal, the excavator was 
fitted with a 1000-L excavating bucket (Fig. 1), which was used 
to spread charcoal and establish the mixing treatments. 

For the charcoal treatments, each plot received 500 + 15 kg 
of biochar, equivalent to 10 tons biochar ha! dry weight, 
which is within the range of addition rates commonly used in 
biochar experiments in agro-ecosystems (Liu et al., 2013). This 
is approximately three times greater than the highest quantities 
of black C that have been estimated in pine forests in northern 
Sweden (Ohlson et al., 2009). The biochar was produced by a 
local company (Vindelkol AB, Vindeln Sweden) and is sold 
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Fig. 1 The excavator used to apply biochar and mix soils in 
an experiment consisting of four treatments; a control, biochar 
only, soil mixing only, and biochar plus soil mixing (n = 6 of 
all treatment combinations). 


and marketed as ‘Terra Preta,’ a product intended as a soil 
amendment (www.vindelkol.se). The biochar was made pri- 
marily from the wood and bark of P. sylvestris, and a small 
portion of Picea abies and Betula pendula. Analysis of the biochar 
showed it had a pH of 8.04 and contained concentrations of 
extractable PO; -P, NO}, and NH; of 1.26, 0.14, and 
1.38 mg ke™!, respectively, and a C concentration of 74%, and 
a P concentration of approximately 300 mg kg~'. Once biochar 
was applied to each plot assigned for biochar addition, the 
excavator was used to mix the soil in all plots assigned for soil 
mixing. This was done with the excavator by removing the soil 
to a depth of approximately 30 cm, piling it within the plot, 
and then evenly spreading and smoothing throughout the plot 
using the steel teeth on the front edge of the excavator bucket. 
The same mixing and spreading procedure was used in the 
mixing only treatment, whereas control plots only received dis- 
turbance associated with tree and stump removal. In the spring 
of 2012, a fence was built around the entire experimental area 
to exclude large herbivores (e.g., moose), and the site was 
planted with approximately 10-cm-tall containerized P. sylves- 
tris seedlings provided by a commercial nursery, with plants 
arranged in a2 x 2 m grid pattern. 


Soil chemistry measurements 


We assessed the relative availability and mineralization of 
NHj , NO; and PO;-P ions by extracting and analyzing bulk 
soil samples and by conducting an in situ mineralization assay. 
In June, 2013, five bulk soil samples were collected from within 
each plot from each of two depths (0-10 and 10-20 cm). These 
samples were collected using a 49-mm-diameter soil corer, 
from the plot center and from the four midway positions 
between each plot corner and the plot center. Samples were 
placed in coolers, returned to the laboratory, and sieved 
(2 mm). The five subsamples from each depth within each plot 
were then bulked, and a 20-g homogenized subsample was 
extracted with 50 mL of 1 m KCI. At the same time these soil 
samples were collected, an additional identical set of 5 cores 
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were collected directly adjacent to the first set, placed in poly- 
ethylene bags, and returned to the same position from where 
they originated. These samples were left in the soil for 
2.5 months, after which they were collected, bulked, processed, 
and extracted in the same way as the original set. Concentra- 
tions of NH] , NO} , and PO; -P were measured on extracts 
of each soil type (i.e., bulk soil and incubated soils) using stan- 
dard colorimetric techniques on an autoanalyzer IN (Omni Pro- 
cess, Solna, SE) (Gundale et al., 2008, 2011a; Wardle et al., 
2013). Approximately 5 g of each soil sample was separated for 
the determination of gravimetric water content, which allowed 
nutrient concentration measures on all soil extracts to be 
reported on a soil dry weight (d.w.) basis. For each plot, net N 
mineralization and P mineralization (mg g7! soil dry weight 
per day) were calculated as the difference in inorganic N (NH; 
plus NO;) or P (PO; -P) concentrations between the nonincu- 
bated and incubated samples, and net nitrification was calcu- 
lated as the difference in NO, concentrations, and divided by 
the incubation period (Eno, 1960; Gundale et al., 2011a). 

In May, 2013, we placed five mixed bed ionic resin capsules 
(PST1 capsule, Unibest, Bozeman, USA) in each plot adjacent to 
the same five sampling locations described above (Gundale 
et al., 2008). The capsules were inserted at a 45° angle, with a 
final depth of 5 cm beneath the surface of the humus, and 
always at least 20 cm from where the soil was sampled. Ionic 
capsules were collected in October 2013 and extracted using 
three consecutive rinsing of 10 mL 1 m KCI (30 mL total) (Gun- 
dale et al., 2011b, 2014). Extracts from these capsules were ana- 
lyzed as described above for NH; , NO}, and PO, -P, and 
nutrient concentrations were reported on a per capsule basis. 

Soil pH was measured on each sample by creating a slurry 
solution consisting of a 1 : 1 mixture of dry soil and deionized 
water. After combining the soil and water, the slurry was 
allowed to equilibrate for 4 h prior to measuring pH. 


Soil microbial community measurements 


We estimated the relative active soil microbial biomass using 
the substrate-induced respiration (SIR) method (Anderson & 
Domsch, 1978), as modified by Wardle (1993) and Gundale 
et al. (2011b). Measurements were made using a composite of 
20 g (d.w.) soil sample from each plot at both depths (i.e., 0-10 
and 10-20 cm). Each sample was placed into a 100-mL glass 
bottle, adjusted to a water content of 125% (d.w. basis), injected 
with a glucose solution, and then sealed. Evolution of CO% 
between 1 and 3 h following glucose addition was then deter- 
mined by injecting 5 mL subsamples of headspace gas into an 
infrared gas analyzer. The difference in respiration between the 
two sample times is an effective proxy of soil microbial bio- 
mass (Anderson & Domsch, 1978). 

In addition to SIR, we used phospholipid fatty acid (PLFA) 
analysis to measure soil microbial community structure (Bligh 
& Dyer, 1959). The PLFAs were extracted from 1 g of freeze- 
dried soil samples using a modification of the Bligh and Dyer 
method (Bligh & Dyer, 1959; McIntosh et al., 2012; Wardle 
et al., 2013). Extracts containing PLFAs were analyzed on a Per- 
kin Elmer Claris 500 Gas Chromatograph (Massachusetts, 
USA), as described elsewhere (McIntosh et al., 2012; Wardle 


et al., 2013). The abundance of identified PLFAs was reported 
as micromoles per gram of organic matter using conventional 
nomenclature and subsequently converted to relative abun- 
dance. Different types of PLFAs represent different compo- 
nents of the soil microflora. The PLFAs 18:1m9 and 18:26 
were used to estimate the contribution of fungi, while the 
branched fatty acids 10me16:0, 10me17:0 and 10me18:0 were 
used to estimate actinobacteria. Bacterial PLFAs including 
i-15:0, x-15:0, 15:0, i-16:0, 16:19, 16:1@7t, i-17:0, cy-17:0, «-17:0, 
18:1m7, cy-19:0 were used to represent total bacteria. Gram- 
positive bacteria were represented by branched fatty acids 
i-15:0, «-15:0, i-16:0, i-17:0 and «-17:0, while cy-17:0, cy-19:0 and 
18:1w7 were used to represent gram-negative bacteria. 


Soil respiration measurements 


We measured total soil respiration (i.e., autotrophic and het- 
erotrophic respiration combined) on six separate measure- 
ment events between May 29 and September 13, 2013, with 
approximately 3 weeks between each event. Each measure- 
ment event required 2 days, with measurements occurring 
between 9:00 and 16:30 hours. The order in which plots were 
measured was changed between sampling events with the 
intention of dispersing any diurnal variation in measurements 
equally among blocks and treatments. The measurements 
were made by installing five cylindrical collars (each 25 cm 
diameter, 10 cm high) in each 0.05-ha plot in early May 2013. 
The collars were inserted 1 cm into the soil surface, and all 
aboveground vegetation inside and 5 cm around each collar 
was manually removed. The collars were allowed to equili- 
brate for 17 days after they were set up before the first mea- 
surement. The height from soil surface to the rim was 
measured in four cardinal directions within each collar to cal- 
culate the headspace volume. Soil respiration (CO, efflux) 
was then measured by sealing the headspace within an opa- 
que plexiglass lid fitted with a portable infrared gas analyzer 
(CARBOCAP model GMP 343, Vaisala, Finland). From each 
collar, CO, concentrations were recorded every 15s for 
3 min. During each measurement, we recorded the headspace 
air temperature using a digital thermometer probe. Soil respi- 
ration within the headspace of each chamber was calculated 
using a linear regression of CO2 concentration vs. time, with 
the slope of the regression indicating the CO, efflux. Esti- 
mated values were subsequently adjusted for variation in 
headspace volume and air temperature, and converted to a 
soil surface area basis, resulting in units of 
umol CO, m~ s™!, as described by Hasselquist et al. (2012). 


Vegetation measurements 


Total vegetation and graminoid cover in each plot were esti- 
mated in August, 2013, by randomly placing 20 0.5 x 0.5 m 
quadrats within the plot and visually estimating percent cover 
values within each quadrat. To assess plant community compo- 
sition within each plot, we used point quadrat analysis in 
August, 2013, as described by Wardle et al. (2003) for each of 
three randomly positioned 0.5 x 0.5 m quadrats. The percent 
cover of all plant species (including vascular plants and 
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mosses) within each quadrat was estimated by downwardly 
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Table 2 The F- and P-values derived from three-factor ANOvAs comparing the main and interactive effects of biochar (added or not 
added), soil mixing (mixed or not mixed), soil depth (0-10 or 10-20 cm) on soil microbial variables, including substrate induced respi- 
ration (SIR), and phospholipid fatty acid (PLFA) data, including total, fungal, and bacterial PLFAs, and the fungal to bacterial ratio. 


Bolded F- and P-values indicate a significant effect at P = 0.05 


Mixing (M)* Biochar (B)* Depth (D)* B x M* B x D* M x D* B x M x D* 

F R Ẹ P F P Ẹ P EF P F P F P 
SIR (ug C g7! h’) 15.2 0.000 25 0.121 129.5 0.000 21.6 0.000 0.6 0.432 17.4 0.000 11.5 0.002 
Total PLFA (nmol g~' O.M.)+ 12.9 0.001 1.7 0.204 83.5 0.000 4.3 0.046 0.2 0.627 169 0.000 2.0 0.170 
Fungi (nmol g~! O.M.) 141 0.001 28 0.102 87.1 0.000 3.3 0.077 0.7 0.399 19.0 0.000 1.6 0.210 
Bacteria (nmol g~! O.M.) 12.0 0.001 1.5 0.226 80.2 0.000 46 0.039 0.3 0.623 15.7 0.000 2.3 0.142 
Gram + (nmol g7' O.M.) 11.6 0.002 1.9 0.183 85.0 0.000 44 0.044 0.3 0.571 15.8 0.000 19 0.174 
Gram - (nmol g`" O.M.) 11.9 0.001 1.3 0.271 73.5 0.000 4.7 0.038 0.2 0.686 15.0 0.000 24 0.134 
Actinobacteria (nmol g~'!O.M.) 12.9 0.001 2.0 0.170 84.0 0.000 4.7 0.038 0.2 0.682 163 0.000 2.0 0.163 
Fungal to bacterial ratio 0.5 0.933 5.2 0.029 38.1 0.000 13 0.266 0.1 0.933 68 0.013 0.1 0.813 


*Degrees of freedom: 1,5. 
+Organic Matter (O.M.). 


Table 3 The degrees of freedom, F- and P-values from a 
repeated-measures ANOVA evaluating the main and interactive 
effects of biochar (added or not added), soil mixing (mixed or 
not mixed), sampling time, and their interactions on soil respi- 
ration, expressed as pmol ms !. Bolded F- and P- values 
indicate a significant effect at P = 0.05 


df F-values P-values 
Mixing 15 7.9 0.006 
Biochar 15 2.8 0.096 
Time 3,15 12.1 0.000 
Mixing*Biochar 1,5 2.5 0.117 
Mixing*Time 3,15 2.8 0.019 
Biochar*Time 15 0.1 0.992 
Biochar*Mixing*Time 3,15 0.8 0.578 


analyzed using SPSS version 21 (IBM, New York, NY, USA) 
except for the PCA for which Primer version 6 software (Luton, 
UK) was used. Data were first evaluated for assumptions of 
normality and homoscedasticity prior to anova and were 
always found to meet these assumptions. 


Results 


Soil response variables 


Numerous soil variables were responsive to biochar 
addition, soil mixing, and soil depth, as well as their 
interactions (Tables 1, 2, and 3). For soil depth, extract- 
able PO; concentrations (Fig. 2a), net N mineralization 
(Fig. 2b), and all soil microbial community parameters 
(Fig. 3) decreased with depth, whereas extractable NO; 
concentrations increased with depth (Fig. 2c). 

Biochar addition caused a significant increase in net 
N mineralization (Fig. 2b, Table 1), extractable NH; 


concentrations (Fig. 2d, Table 1), resin-sorbed NHj con- 
centrations (Fig. 4a, Table 1), and the fungal:bacterial 
ratio (Fig. 3e, Table 2). Biochar had no significant main 
effects on any other response variable nor were there 
any significant interactions between biochar and soil 
depth. 

Soil mixing had a significant main effect on many 
response variables. Mixing caused a significant 
decrease in net N mineralization rates (Fig. 2b, 
Table 1), SIR, total PLFAs, the abundance of all 
microbial functional groups (Fig. 3a-d, Table 2), soil 
respiration rates (Supplementary Fig. 1, Table 3), and 
an increase in soil NO; concentrations (Fig. 2c, 
Table 1). For several variables, soil mixing also 
showed an interactive effect with soil depth (Table 1, 
2). For net N mineralization, this interaction occurred 
because the negative effect that soil mixing had on N 
mineralization was more pronounced in the surface 
compared to deeper soil (Fig. 2b). Likewise, for all 
microbial response variables (Table 2), the negative 
impact of mixing was more pronounced in surface 
relative to deeper soil (Fig. 3a-e; data for Gram+, 
Gram-, and actinobacteria not graphically depicted). 
For soil pH, the interactive effect occurred because 
mixing increased soil pH in the surface 
(pH = 4.88 + 0.054 and 5.17 + 0.036 in the unmixed 
and mixed soil, respectively), but not in the deeper 
soil layer (5.41 + 0.036 regardless of treatment). For 
soil respiration, the effect of mixing interacted signifi- 
cantly with time (Table 3), with negative effects more 
pronounced in the middle compared to the beginning 
or end of the growing season (Fig. S1). 

Numerous soil variables responded to the interactive 
effect of biochar and mixing, and one additional vari- 
able (SIR) responded to a three-way interactive effects 
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Fig. 2 The mean (+SE) extractable soil PO; concentrations in response to soil depth (a), net N mineralization rates in response to 
soil depth and biochar and soil mixing treatments (b), NO; concentrations in response to soil depth and soil mixing (c), and soil NH 
in response to soil depth and biochar application (d). Within each panel, different capital letters (A or B) above bars or bar groups 
indicate a significant difference between the soil depths; lower case letters (a or b) indicate significant differences across all bars as 
indicated by post hoc Student-Neuman-Keuls tests, and stars (*) indicate significant differences between the two bars within each soil 


depth, as determined using post hoc Student t-tests. 


of biochar, mixing, and soil depth (Table 1 and 2). For 
net N mineralization (Fig. 2b), all microbial parameters 
except fungi and the fungal to bacteria ratio (Fig. 3a-e), 
and resin-sorbed NO, and H2PO, (Fig. 4b,c), two-way 
interactive effects between biochar and mixing occurred 
because the impacts of soil mixing on these variables 
were less severe when biochar was also added. Like- 
wise, the three-way interactive effect found for SIR was 
due to the negative effects of soil mixing being less 
severe when biochar was also added, but only in the 
surface soil (Fig. 3a). 


Vegetation response variables 


The addition of biochar had a significant main effect 
on only one vegetation variable, that is, plant species 
richness, which declined when biochar was added 
(Fig. 5d, Table 4). Soil mixing had a negative impact 
on total plant cover and graminoid cover (Fig. 5a, 
Table 4). Mixing also reduced PCA axis 2 scores 
(Fig. 5b, Table 4), where lower scores represented 
greater abundance of Vaccinium vitis-idaea and Calluna 


vulgaris, and higher scores represented greater 


abundance of Deschampsia flexuosa and Epilobium angus- 
tifolium (loading scores —0.44, —0.22, 0.34, and 0.76, 
respectively). In contrast, pine seedling survival was 
significantly enhanced by soil mixing (Fig. 5c, Table 4), 
as was the germination percentage of all four tree spe- 
cies (Fig. 6, Table 5). For tree seed germination, the 
positive effect of mixing on germination differed 
among the four species (i.e., indicated by a significant 
mixing by species interaction; Table 5). These interac- 
tive effects occurred because the positive effect of mix- 
ing was less pronounced for B. pendula compared with 
the three conifer species (Fig. 6). One variable, that is 
total vegetation cover, was also affected by the interac- 
tion of biochar and mixing, whereby the negative 
effect that mixing had on total cover was more severe 
when biochar was added compared to when it was 
not added (Fig. 5a). 


Discussion 


Consistent with our first hypothesis (i.e., that addition 
of biochar would enhance nutrient availability and min- 
eralization rates, microbial biomass, and soil respiration 
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Fig. 3 The mean (+SE) substrate-induced respiration (SIR) rate (a), concentrations of total (b), fungal (c), and bacterial (d) phospho- 
lipid fatty acids (PLFAs), and the fungal:bacterial ratio (e) in response to soil depth (0-10 or 10-20 cm) and a factorial combination of 
biochar application and soil mixing. Within each panel, different capital letters (A or B) above each group of bars indicates a signifi- 
cant difference between the soil depths, whereas lower case letters (a or b) indicate significant differences across all bars, as indicated 


by post hoc Student-Neuman-Keulls tests. 


rates independently of the soil mixing treatment), we 
found that biochar increased NHj availability, regard- 
less of whether it was applied to the surface or mixed 
into the soil. However, a majority of other variables 
were unresponsive to biochar addition, such as NO; 
and PO,-P, most microbial community variables, and 
soil respiration. One potential explanation for the 
observed increase in NHj concentrations is that the bio- 
char itself served as a source of NHj , which accumu- 
lates within ash residues on its surfaces during the 
charring process (Gundale & Deluca, 2006; Pluchon 


et al., 2014). However, if we consider that the mean soil 
bulk density at the site was 1 g cm © and a depth of 
20 cm, the increase in NH; concentrations caused by 
biochar addition scales up to approximately 0.5 kg NH} 
per plot, which is 724 times greater than the NH} 
directly contributed by the biochar itself (i.e., 500 kg of 
biochar per plot x 138mg NH} kg! bio- 
char = 690 mg NH; added per plot). Thus, the positive 
effect of biochar on soil NHj concentrations was instead 
likely to have arisen by promoting net N mineralization 
rather than ash input (indicated by a main biochar 
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Fig. 4 The mean (+SE) resin-sorbed NH; in response to bio- 
char addition (a), and resin-sorbed NO, (b) and PO, (c) in 
response to biochar application and soil mixing treatments. 
Within each subpanel different lower case letters above bars, or 
stars (*) above bar pairs indicate significant differences, as 
determined using post hoc Student's t-tests. 


effect, Table 1), which was particularly evident in the 
mixed plots (indicated by a biochar by mixing interac- 
tive effect, Table 1). Enhancement of N mineralization 
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associated with biochar especially in the mixed plots 
may have occurred as a result of enhanced gas 
exchange or soil moisture retention associated with its 
very high micropore volume (Keech et al., 2005; Pluchon 
et al., 2014). 

Contrary to the first part of our second hypothesis 
(i.e., that soil mixing alone would enhance nutrient 
availability and mineralization rates, and soil respira- 
tion, but decrease total microbial biomass and fungal 
to bacteria ratios), we found that soil mixing caused a 
decline in a majority of these variables. While all plots 
within our experiment were subjected to a relatively 
minor level of disturbance during forest harvest and 
stump removal (which left the organic horizon mostly 
intact), our mixing treatment caused a substantially 
more severe disturbance by turning over the entire 
soil volume to a depth of approximately 30 cm, 
thereby likely affecting a wide range of soil properties 
such as aggregate structure, porosity, and gas 
exchange. The disturbance caused by the mixing treat- 
ment resembled several mechanical site preparation 
techniques (e.g, mounding or furrowing) widely 
employed in the boreal region used to enhance forest 
establishment, which have been shown in numerous 
studies to enhance short term nutrient availability and 
mineralization rates, as well as bacterial biomass, and 
soil respiration rates (Frey et al., 2003; Siira-Pietikainen 
et al., 2003; MacKenzie et al., 2005; Piirainen et al., 
2007). These changes can occur because soil distur- 
bance causes soil organic matter stabilized within soil 
aggregates to become more available to the microbial 
community, leading to short-term increase in bacterial 
biomass, soil respiration, nutrient transformation rates, 
and nutrient accumulation in the soil (Siira-Pietikainen 
et al., 2003), followed by an eventual decline (Piirainen 
et al., 2007). One potential explanation for the contra- 
dictory results found in our study is that these effects 
are short lived and had ended prior to our sampling 
effort. 

We found no evidence for our second hypothesis that 
soil mixing would enhance effects of biochar on mea- 
sured soil variables. Instead, we found that biochar 
dampened many of the impacts that mixing by itself 
had on soil nutrient availability, nutrient transformation 
rates, and microbial community variables. Several prop- 
erties of biochar may have contributed to this dampen- 
ing. For instance, numerous studies have shown that 
biochar can enhance microbial biomass and activity by 
serving as a substrate and habitat for microbes, and by 
improving soil physical properties (Pietikainen et al., 
2000; Lehmann et al., 2011). Our finding of enhancement 
of SIR and soil bacterial and fungal PLFAs in the mixed 
soil with vs. without biochar would have likely led to 
the enhancement of net N mineralization rates, given 
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Fig. 5 The mean (+SE) graminoid and total cover of the ground layer vegetation in response to soil mixing and biochar treatments 
(a), the mean (95% confidence interval) of the first and second axis of a principle component analysis (PCA) describing the vegeta- 
tion composition (b), and mean (+SE) survival of P. sylvestris seedlings (c) in response to soil mixing treatments, and the response on 
species richness to the charcoal treatment (d). Different letters above bars (a or b) in subpanels a, c, and d indicate significant differ- 


ences, as determined using Student-Neuman-Keuls post hoc tests (a) or analysis of variance (c and d). 


Table 4 The F- and P-values from anova evaluating the main 
and interactive effects of biochar (added or not added) and soil 
mixing (mixed or not mixed) on plant species richness, PCA 
axis 1 and 2, graminoid and total cover, and pine survival. Bol- 
ded F- and P- values indicate a significant effect at P = 0.05 


Mixing (M)* Biochar (B)* B x M* 

Ẹ P F P FP 
Species Richness 2.1 0.166 85 0.011 0.1 0.775 
Total cover (%) 62.4 0.000 3.3 0.088 5.5 0.033 
Graminoid cover (%) 117.3 0.000 0.8 0.394 1.7 0.218 
PCA 1 1.0 0.345 0.5 0.495 0.0 0.960 
PCA 2 84 0.011 0.0 0.971 1.1 0.314 
Pine Survival (%) 145.9 0.000 0.2 0.630 2.5 0.137 


*Degrees of freedom: 1,5. 


this process is driven by the soil microbial community. 
This stimulation of the soil microbial community may 
also have immobilized NO}, thus leading to the lower 
NO, concentrations we observed in the mixed soil with 
vs. without biochar. The significant increase in soil 
PO;-P in mixed soils with vs. without biochar could 
have been due to biochar interfering with PO, -P com- 
plexation into humic substances, or with Fe and Al ions 
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Fig. 6 The mean (+SE) germination rate (%) for planted seeds 
of Betula pendula, Picea abies, Pinus sylvestris, and Pinus contorta 
in response to soil mixing 92 days after seed sowing. Stars 
above each pair of bars indicate a significant difference, as 
determined using post hoc Student's t-tests whereas, capital let- 
ters above each bar pair indicate a significant difference in ger- 
mination rate across the four species, as determined using 
Student-Neuman-Keuls post hoc tests. 


(Lehmann et al., 2003; Topoliantz et al., 2005), which can 
limit P availability in the boreal region (Giesler et al., 
2005). 
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Table 5 The degrees of freedom, F- and P-values from a split 
plot anova evaluating the main and interactive effects of bio- 
char (added or not added), soil mixing (mixed or not mixed), 
seedling species, and their interactions on the percent germina- 
tion of four seedling species (Betula pendula, Picea abies, Pinus 
sylvestric, and Pinus contorta). Bolded F- and P-values indicate a 
significant effect at P = 0.05 


df F-value P-value 
Mixing 1,20 23.4 <0.001 
Biochar 1,20 2.4 0.131 
Species 3,441 187.6 <0.001 
Mixing*Biochar 1,20 0.109 0.745 
Biochar*Species 3,441 0.070 0.976 
Mixing*Species 3,441 5.215 0.002 
Mixing*Biochar*Species 3,441 0.403 0.751 


Contrary to our third hypothesis (i.e., that treatment 
effects on nutrient availability would correspond with 
differences in plant community composition), we found 
that elevated soil NH{ concentrations that occurred in 
response to biochar addition did not correspond with 
large changes to the ground layer vegetation. The only 
effect of biochar was to cause a small decrease in plant 
species richness per unit area. However, in partial sup- 
port of our third hypothesis, we found that mixing 
reduced graminoid and total plant cover and promoted 
dominance by ericaceous shrubs; this corresponded 
with the lower net N mineralization rates and resin- 
sorbed PO, -P concentrations in the mixed plots. The 
most likely explanation for the observed unresponsive- 
ness of vegetation to biochar application, as well as for 
the relatively weak relationship between vegetation 
cover and nutrient availability, is that vegetation 
appeared to be primarily controlled by disturbance 
rather than soil nutrient availability over the time scale 
of our study. The soil mixing treatments clearly dam- 
aged and buried vegetation, resulting in a direct reduc- 
tion in vegetation cover. While no effect of enhanced 
nutrient availability on vegetation was evident in our 
study, we anticipate that these effects could emerge 
over longer times scales, as the relative importance of 
plant competition for nutrients gradually becomes more 
important (Grime, 1979). 

Contrary to the predictions of our fourth hypothesis, 
our data showed that biochar application had no dis- 
cernable short-term impact on tree seed and seedling 
response despite significantly increasing extractable 
NHj concentrations. This result is consistent with some 
previous studies in northern conifer forests (Naydenov 
et al., 2006; Heiskanen et al., 2013) but not others (War- 
dle et al., 1998; Nilsson et al., 2008; Pluchon et al., 2014). 
Pluchon et al. (2014) showed that biochar effects on 
seedling growth of boreal forest tree species were 
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determined by interactive effects between biochar and 
soil, whereby biochar with higher concentrations of 
available PO, -P had large effects only in soils where P 
was limiting relative to N. This result could explain the 
unresponsiveness of seedlings to biochar application, as 
biochar may not have alleviated the primary resource 
limitation or reduced competition with the ground layer 
vegetation at the site over the time scale of our study. 
However, in support of our fourth hypothesis, P. sylves- 
tris seedling growth and germination rates of four sown 
species (P. sylvestris, P. contorta, P. abies, and B. pendula) 
were strongly impaired in the plots where vegetation 
cover was the greatest (i.e., the unmixed plots). These 
results highlight the importance of ground layer vegeta- 
tion in determining tree establishment success in boreal 
forests (Thiffault et al., 2013). 

The suggestion that biochar management can be used 
to effectively mitigate global anthropogenic CO, emis- 
sions require its successful implementation in nearly 
every major biome, and thus, far studies evaluating bio- 
char application in forest environments are almost com- 
pletely lacking. Our study provides several insights into 
understanding the impacts and effectiveness of biochar 
management in forested environments and boreal for- 
ests specifically. First, it suggests that biochar applica- 
tion can enhance availability of some soil nutrients 
known to limit productivity in the boreal region (nota- 
bly NHj;) while directly contributing to a relatively sta- 
ble pool of soil C. Given that we applied 10 t ha! of 
biochar with a C content of 74% (ie., 7.4 Mg C ha!) to 
a soil containing approximately five times as much C 
(i.e., 38.8 Mg ha~’), and no changes in soil C efflux (i.e., 
soil respiration) were detected, our study suggests that 
biochar was relatively stable and at most caused only 
minor losses of existing soil C. This finding suggests 
that biochar may indeed serve as a useful tool to simul- 
taneously enhance soil fertility and soil C stocks in bor- 
eal forests. Secondly, our data provide new insights into 
the potential impacts of biochar on aboveground C 
sequestration in the boreal region. Biochar alone was 
ineffective at enhancing growth of the ground layer veg- 
etation or the regeneration of forest trees over the 2-year 
duration of our study, contrary to a wide range of stud- 
ies conducted in agro-ecosystems that have shown posi- 
tive effects on plant growth (Jeffery et al., 2011). Given 
that trees eventually account for most of the above- 
ground C in establishing forest stands, our results sug- 
gest that biochar application in the boreal region might 
be more useful and successful as a C storage and 
sequestration tool if used in combination with site prep- 
aration systems aimed at optimizing forest regeneration. 
Within this context, biochar may help prolong the short- 
term transient effects that these preparation techniques 
may have on soil nutrient availability, microbial 
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biomass, and activity, while directly enhancing recalci- 
trant soil C stocks. While our study covered a relatively 
short response period relative to a typical forest rotation 
length, it is a critical first step in evaluating the impacts 
of biochar management in forested ecosystems on the 
global C cycle. Our study, in combination with addi- 
tional long-term studies, is necessary before biochar 
management can be accepted and promoted within pro- 
posed C trading schemes in the boreal region. 
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Figure S1. The mean (+SE) soil respiration rates in 
response to soil mixing, measured at six sampling times 
during the growing season, 2013. 
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